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Abstract 

The Friedmann-Robertson- Walker metric with spherical topology is 
(**] calculated as an effective metric at the brane in a multiple branes in 

D— dimensional spacetime scenario. In this model the radius of the brane 
is the cosmological scale factor, and its evolution is calculated from the 
y—i viewpoint of the observers on their respective branes. The cosmology of 

the brane-universe is analyzed in cases where the anisotropic pressure has 
a cosmological constant like state equation. The pressure values needed 
to distinguish the oscillating, expanding or collapsing solutions are found, 
l^"* and the minimum value of a brane-universe mass to prevent a collapse into 

,* ^ a black hole is also calculated. Finally the equation of motion is solved 

(3jT) numerically in order to illustrate the different cosmological scenarios and 

i i to validate the analytical results. 
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1 Introduction 



The Fricdmann- Robertson- Walker metric (FRW) describes a homogeneous isotropic 
universe that can be modeled as a perfect fluid. The dynamics of this fluid is 
given by a cosmological scale factor, which determines the spatial volume of the 
universe, and depends on the amount of matter in the universe [l]. This model 
became the standard cosmological model due to the observations of Hubble in 
1929 1 2] showing an expanding universe in opposition to predictions of the Ein- 
stein's static universe. This observation led to the conclusion that the universe 
was hotter and denser in the past with an infinite density at the Big Bang. Al- 
though such initial singularity causes difficulties to a rational description of the 
universe, that model gained prominence in the 70's, due to the Hawking and 
Penrose's singularity theorems (31. However, the latest observations indicate 
that the universe is expanding at an increasing rate (4j[5] what requires some 
kind of matter which gravitate repulsively like a cosmological constant. This 
kind of matter violates the conditions of the singularity theorems, making its 
results inapplicable to the nowadays observed universe. Since then, models with- 
out cosmological singularity returned to the scene as possible candidates for the 
description of the universe like for example models of oscillating universe (6j|9| . 

Brane-uni verse models emerged to solve theoretical problems, like the uni- 
fication of fundamental interactions and the hierarchy problem (10-15 . These 



models describe the universe as a flat brane, and because of this it can not 
be used to solve cosmological problems, like the accelerated expansion and the 
initial singularity. Thus, other geometries are necessary but only two are com- 
patible with the FRW metric: the spherical and hyperbolic ones. For easily 
incorporating the observed symmetries of the universe, like the galaxy isotropic 
runaway, the existence of a preferred frame, and a cosmic Gaussian time, the 
spherical brane model has received great attention. The phenomenology of 
these models has been studied in past decade (16 17 and for more recent stud- 
These models has been proved compatible with the observational 
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data 21-23 , as well as, by introduction of different cosmological constant in 
the bulk, with the observed dynamics, modeling the dark energy |24| . 

In this piece of work we continue a recent study where a multiple anisotropic 
concentric spherical branes model was constructed |25 |. Here, a dimensional 
reduction is made to obtain the effective metric measured by observers on the 
brane and the evolution of the cosmological scale factor in the description of 
these observers. This reduction give us the cosmology of the brane-uni verses, in 
such a way as to be compared with the current data. A more detailed study of 
the cosmology is made for cases where the anisotropic pressure is of cosmological 
constant type. This becomes important since this model naturally tends to 
expand, in agreement with the current observations. Studying the null velocity 
points we can found the necessary conditions ( but not sufficient) to models 
of oscillating or eternal expanding universes. In the oscillating scenario we 
analyze the conditions able to prevent the collapse of the universe into a black 
hole. To reinforce the study of null velocity points the equation of motion is 
solved numerically. At this point we will discuss the energy condition in the 
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oscillating scenario. 

This work is organized as follows: Sec. II gives a short review of results ob- 
tained in multiple spherical brane model. In particular the metric, the equation 
of motion and the possible state equation for anisotropic pressure are shown 
since they are needed for the sequence of this work. In Sec. Ill the effective 
metric on the brane is calculated through a dimensional reduction, and the cos- 
mology of the brane to observers on them is also obtained. In Sec. IV the 
possible models of a universe with state equation for anisotropic pressure like 
a cosmological constant is discussed, and the conditions for the existence of 
oscillating, collapsing or eternally expanding universes are studied. In Sec. V 
the equation of motion is solved numerically to illustrate the behavior of the 
cosmological scale factor in different scenarios, as well as to confirm the results 
obtained in Sec. IV. Also in this section we discuss the weak energy condi- 
tion in oscillating scenario. The last section shows the conclusion and possible 
perspectives of this work. 



2 The spherical brane model 

Following the multiple spherical brane model developed in ref. [25] the Einstein's 
equation is solved in a D dimensional bulk with spacetime dependent cosmolog- 
ical constant. The matter field has n spherical and anisotropic (D — 2)-branes, 
and the obtained invariant line element is 

dr 2 = -A(r, t)dt 2 + B(r, t)dr 2 + r 2 dn 2 D _ 2 , (1) 

where Qd-2 is the angular line element in D dimensions, formed by D — 2 
angular variables. The metric elements are 

B -, (r:t) = 1 _^Mpl_ r ^ t) (2) 
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where Pi is the anisotropic pressure, 2AA^ = (D — 1)(D — 2)AAi is the difference 
between the internal and external cosmological constant of the z-th brane, and 
for completeness Mq = Rq = 0. Since the branes are dynamic, the following 
equation of motion for the z-th brane is obtained from the conservation of energy- 
momentum tensor: 



' dU 



-Bi(Ri) [Pi + Pi -2 Pi Uf] 
, G D M{R % ) 



[Pi - Pr (7i + U?)] 
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R 
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where dti = A i (R i ) / B i (R i )dt, L/j = and 7$ is the ratio between the 
pressure and energy density of matter on the brane. The achievement of the 
above equation of motion gives the following two state equations that remove 
divergences 

Pi = -Pi or Pi = Pi Uf. (9) 



With the above anisotropic pressures the brane can have any value of matter 
on it. This freedom allows us to keep the cosmological eras, in same way the 
standard cosmological model. 



3 Effective metric and cosmology on the brane 

In order to calculate the effective metric and its evolution measured by observers 
on a specific brane a reduction is performed to describe the brane-universes 
cosmology. To do this we will start with the D-dimensional metric ([lj, and 
using the relationship r — Ri (t) . we can make the dimensional reduction on the 
brane to obtain the effective metric 



dsi = - [MR,) - Bi{Ri)V^] dt 2 + R t (t) 2 dn 2 D _ 2 



(10) 



To simplify the above expression, we use the following comoving time and radial 
coordinate on the brane 



dr, 



Ai{Ri)- Bi{Ri)V 2 dt, 
r = sin 9q . 



(11) 
(12) 
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In this coordinates, the reduced metric takes the form 

" df 2 



dsj 



1-f 



f 2 dn 2 D _ 3 



(13) 



This is the Friedmann-Robertson- Walker metric for a universe with (D — 2) 
spatial dimensions with spherical topology, where the brane radius is the cos- 
mological scale factor. To obtain the scale factor evolution, Ri(ri), we need 
to calculate the evolution law in terms of the new variable r,. Using the 
composed derivation rules we obtain 



Ui 
dXh 



B ^ Wi 

1 + Bi{IU)Wf " 



Bi(Ri) 



(l + B^W 2 ) 2 
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(14) 
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where Wi = dRi/dri and the prime indicates the derivation with respect to the 
argument. In terms of comoving time the equation of motion (pi takes the form 



Pi 



dm 
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b {Ri) + b 2 {R l )w 2 + h(Ri)wt, 



(16) 



where, to simplify, we defined w; = J Bi(Ri)Wi. The above coefficients for the 
P L = —pi case are 
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and, for Pi = PiU 2 , the coefficients are 



Bi(Ri) 
G D M{Ri) 



R 



D-2 



Bi(Ri 
XiRi , 



A,; Ri 



(17) 



bo(Ri) 
b 2 {Ri) 



AA, 



K D Bi(Ri) 



(D-2) 



Bi(Ri)Ri 



{D-2,) 



G D M(Rj 
R?- 2 



— RiK 



AA, 



| 2(D-2) llPl 

Bi(Ri)Ko Bi(Ri)Ri 
{D - 2)pcu 
Bi(Ri)Ri 



(18) 



4 



Equation ( 16 ) gives us the brane-uni verse cosmology as felt by observers on the 
brane, and as one can see, from velocity independent coefficient, the cosmologies 
are strongly influenced by the choice of radial pressure. In the first case (Eqs. 
(17)), the brane-uni verse tends to expand due to the negative pressure ( like 



a cosmological constant), while the in the second case (Eqs. 
collapse. 



(18)) it tends to 



4 Cosmology for Pi = —pi 



In this section we study the model with a cosmological constant like state equa- 
tion. We choose the P, : = — pi case because it has a state equation that simplifies 
the analytical calculations. Furthermore, this state equation is consistent with 
the energy-momentum tensor used in the model, and provides a natural explana- 
tion for the condition expressed in Eq. ^ . There are two different cosmological 
scenarios depending on the signal of acceleration when the velocity of the i-th 
brane vanishes. That happens for Wi — 0, and the equation (16 1 give us the 
acceleration 



AA,; 



w,=o 



KDBi(Ri) 



(D - 2) pi 



l+7j 
Bi(Ri)Ri 



(19) 



Considering only the region where the energy density, pi, is positive, the accel- 
eration will be positive if 

AA, D - 2 , 
< - + —5—{^ + li)pi 

K D Hi 

oc AXiRP - 1 [(£> - 1) + (D — 2) (1 + 7i )] + 

+2{D-2){l + ll )G D M l . (20) 

Assuming that 7,; > — 1, what is physically consistent since we don't know 
anything that violates this condition, the acceleration in the rest points will be 
positive if 

2(D-2)(l+ li )G D M i D _, 
-p-l) + P-2)(l + 7l ) <AA ^ ' (21) 
For the AAi > case, the solution is trivial (Ri > 0), and this universe model 
will expand indefinitely, which makes it consistent with the observations. For 
the AA; < case we obtain the non-trivial solution 



R, < 



2G D Mj 
IAAJ 



D - 1 

(D-2)(l +ll ) 



-1 l/D-l 



(22) 



The above result shows that, in rest points, the acceleration will be positive 
if the brane radius is smaller than the above quantity, and vice-versa. The 
obtained result defines a critical radius 



R'; 



2G D Mj 
lAXAai 



1/(23-1) 



(23) 
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where 

The adopted procedure don't give us the return points, but show the existence 
of a rest point with negative acceleration beyond the critical radius and other 
with positive acceleration before it. This condition is required for an oscillating 



universe model. However, the equation of motion ( 16 ) is valid just beyond the 
Kottler horizon. Therefore, in order to have an oscillating universe is necessary 
that the critical radius be bigger than the horizon radius, 



2G n Mi 



l/C-D-l) 



* < [wS) ' (25) 

Since the horizon is defined by the pole of the metric (J2|, given by equation 

- \i (R^f- 1 + (Rtf- 3 2G D ]T Mj = 0, (26) 

3=0 

we can find the mass that makes the critical radius outer the horizon, which is 

-(lAA.I^ + A.OIAA.K-i-I, (27) 

where 



-l 

M™\ = Y. M 3- ( 28 ) 

3=0 



The above expression has analytical solution just for the most inner brane, where 



M™\ = 0. To this brane the critical radius will be bigger than the horizon only 
if 



Mi < J^^cti (|AAi|<7i + Ai)~ (£, ~ 1)/2 . (29) 

This is the maximum mass where the most inner brane can have in order to not 
collapse. If the mass is bigger than this value the evolution law forces that it is 
contracting and pass through the horizon creating a black hole. If the mass is 
less than its maximum value the brane can oscillate as in bouncing models. 



5 Numerical Results 



In this section the equation of motion (Eq. (16l) is solved numerically. The 
parameters values used in the calculations were chosen in a way to show the 
general behavior expected according for different scenarios shown in last sec- 
tion. We should point that the numerical results are more complete than the 
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Figure 1: Numerical solution for D = 5, G5M = 2 and K5 = 1. The dashed 
line is the solution for dust matter on the brane while the solid line is for a 
photon gas. Here we considered an expansive scenario with Ao = Aj = and 
i?i(0) = 100 and Wj(0) = as initial conditions. 



analytical ones. In this case we have found just a necessary condition to dis- 
tinguish between the different scenarios, but here we find the complete solution 
once the initial conditions are given. 

Fig. [T] shows the result for a scenario without cosmological constant, i.e., 
an expanding universe model. The curve shows the behavior of the radius R 
against time. From the figure, one can see an expansionary phase preceded by 
a collapsing stage indicating a model with bouncing. For the used parameters 
the universe expands rapidly with time. In Fig. [2] we consider a scenario of 
an oscillating universe, in accordance with the parameters obtained previously. 
The result shows that the oscillation amplitude must be limited. This hap- 
pens because when considering a negative difference of cosmological constants 
(needed for an oscillating scenario) the energy density on the brane becomes 
non-positive defined. This can be seen from the definition Q 



-^—-^2 \?G D Mi - |AAi|i?f _1 ] 



(30) 



In order to ensure the weak energy condition 26 we have to limit the brane 
radius to 



R i. < 



2GpMj 
IAAJ 



l/D—l 



(31) 



so, we can overestimate the oscillation amplitude 
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Figure 2: Numerical solution for D — 5, G$M — 2 and k& = 1. The dashed 
lines are the solutions for dust matter on the brane while the solid line is for 
a photon gas. we simulate an oscillating scenario, for this we use Ai = and 
A = 4 x 1(T 4 whit Ri(0) = 9.0 and Wi(0) = as initial conditions. 




t 



Figure 3: In above simulation we use the same values of Fig. [2] but we fix 
i?i(0) = 9.99 and 7 = 1/3. The dashed line indicates the critical radius R^. 
This illustrate the different behavior for the regions below and above the critical 
radius. 



assuming the strong energy condition, 7.; < 1/3, the above result leads us, in 
D = 5, to 

(33) 

The parameters used in the simulation give us the value A, < 3.2. Fig. |] 
shows an amplitude of w 0.8, within the estimated range. The difference occurs 
because we considered Ri(0) = 9.0 ( while pi — > when Ri — > 10.0) and the 
behavior for the region Ri > Rf is different of the behavior for Ri < R^, as we 
can see in Fig. [3] where results for Ri(0) — 9.99 provides us with an amplitude 
of 2.0, which is still much lower than the overestimated result. These results 
show the validity of the analytical study of the previous section. 



6 Conclusions and Perspectives 

In this work we compute the effective metric measured by observers on the 
brane. We obtained, in agreement with the standard cosmological model, the 
Friedmann-Robertson-Walker metric with spherical topology. We also obtained 
the cosmology of these brane-uni verse, that describes the evolution of the cos- 
mological scale factor. We observed that the brane cosmology is affected by 
bulk cosmological constant, by matter on the brane, as well as by the brane 
state equation. The matter dependence give us the idea of cosmological eras, 
which depends on the state equation of this matter. 

We also study the cosmology obtained from evolution equation ^ in the case 
where the anisotropic pressure has a cosmological constant like state equation. 
We show that this kind of brane tends to expand, due the negative pressure, so 
it is able to describe the universe without the introduction of cosmological con- 
stants in the bulk. In case where the difference of cosmological constants in the 
bulk is positive or null we found an eternal expanding universe. An interesting 
scenario was found using an negative difference between cosmological constants. 
This has the effect of inhibit the expansion and we obtained an oscillating uni- 
verse. In this scenario we also obtain the mass limit which a brane could have 
in order to not collapse into a black hole. Finally, we solve numerically the 
equation of motion in order to illustrate and reinforce the analytical results. In 
the oscillating case we estimate the oscillating amplitude imposing the strong 
energy condition in the brane. 

An extension of this work is the study of cosmology generated by the velocity 
dependent state equation and study the phenomenology of the model. Despite 
its advantages, we should point that the spherical brane scenario share some 
problems of plane brane models. One of this problems is related to the localiz- 
ability of fields on the brane (27-32 . This can have important phenomenologi- 



cal implications to the standard model and cosmology |33f|35| . However, in the 
present manuscript we have supposed that the matter is localized in the brane 
and we leave the study of trapping of fields in the spherical brane scenario for 
a future work. 
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